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1 Introduction 



The decays K L — > 7r°77, K L — > 7r°e + e~7 and K L — > 7r°7r°7 can all be used as tests of 
chiral perturbation theory. In particular, predictions for the branching ratios of these 
modes show significant increases when one uses 0(p e ) versus 0(p 4 ) chiral perturbation 
theory. The first measurements of Kl — > 7r°77 p] |2j were a factor of three higher 
than the 0(p A ) prediction, but were consistent with the 0(p 6 ) calculation. [3] This 
was also seen in the decay Kl — > 7r°e + e~7 decay where the 0(p A ) prediction was 
inconsistent with the measurement, but consistent with the 0(p 6 ) calculation, jl] For 
the Kl — > vr°7r°7 decay, the branching ratio vanishes at 0(p A ) chiral perturbation 
theory, yet is non-zero at higher order. Recent predictions for this decay range from 
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In addition, the two decays Kl — > vr°77 and Kl —> 7r°e + e _ 7 are important 
for understanding the direct CP violating decay Kl —> 7T°e + e~. Three components 
contribute to the Kl — ► 7r°e + e~ amplitude: direct CP violation, indirect CP violation 
(and an interference term), and a CP conserving term. Recent measurements of the 
decay Kg — > 7r°e + e~|Zj and K$ — > vr / u + /i~jH] have helped to determine the indirect 
CP violating contributions to Kl — > 7r°e + e~ and Kl — > vr°/i + yU _ . The magnitude of 
the CP conserving contributions to Kl — > 7T°/I can be determined by measurements 
of the decay K L —> 7r°77jHJ [TD| and K L —> 7r°e + e~7. The CP conserving term is 
estimated to be small [TT]. 



2 The KTeV Experiment 

The KTeV experiment, shown in Figure Q is a Fermilab fixed target experiment. The 
detector ran in two different configurations: E799 and E832. The E799 configuration 
took advantage of a higher kaon flux to search for rare kaon decays. The E832 
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Figure 1: The KTeV detector. 



configuration was used primarily for a measurement of e'/e. The main difference 
between the two configurations was the use of a regenerator to produce K$ decays 
in the E832 configuration. The two experimental configurations also had a few other 
differences as noted below. 

The KTeV detector contains a charged spectrometer with four drift chambers, 
two on either side of a large dipole magnet. At the downstream end of the detector 
is a two-meter square calorimeter consisting of 3100 pure Csl blocks. Following the 
calorimeter are 10 cm of lead and 5 meters of steel which act as a muon filter. Two 
planes of scintillator, used for muon detection, are located just downstream of the 
steel. Photon vetoes to detect the the presence of particles that would otherwise 
escape detection surround the spectrometer. The charged spectrometer achieves a 
hit resolution of better than 100 /an, while the Csl calorimeter obtains better than 
1% energy resolution over the range of energies of interest. Just upstream of the 
Csl calorimeter is a transition radiation detector (TRD) capable of e/ir separation 
of 200:1 with a 90% efficiency. The TRDs were employed during E799 running, and 
were moved out of the way during E832 running. 

The KTeV experiment took data during a number of different periods between 
1996 and 1999. In the E832 configuration, we had three different running periods: 
1996, 1997 and 1999. The E799 running occurred in 1997 and 1999. Between the 
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1997 and 1999 runs, a number of upgrades were made to the detector to increase 
its reliability and to improve its live time. In addition during E799 running, the 
transverse kick from the magnet was reduced from 205 MeV/ c to 150 MeV/c, enabling 
a larger acceptance for high multiplicity decay modes. For the entire E799 data set, 
approximately 6.6 x 10 11 kaon decays occurred in the KTeV detector. This large kaon 
flux allows us to have an unprecedented sensitivity to a number of rare kaon decays 
with large multiplicity final states. 
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Figure 2: The kaon reconstructed decay position. The green histogram is the Kl 
7r°77 Monte Carlo, while the black histogram is the Kl — > 7r°7r 7r° Monte Carlo. 
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3 Measurement of Kl —> 7r 77 

In the decay Kl — > 7r°77, the final state consists of four photons, two from the 7r° 
and two from the Kl- Our analysis requires four electromagnetic clusters in the Csl 
calorimeter, with no tracks or extra clusters. One can combine the four clusters in 
three different combinations. We choose the combination in which the two photon 
invariant mass reconstructs closest to the known 7r° mass. 

Backgrounds to this decay originate from two main sources, Kl — > 7r°7r° and 
Kl — > 7r 7r°7r°. The first background has the same topology as our signal events. 
However, the relatively small branching ratio for Kl —* 7r°7r° and kinematic cuts 
reduce this background to a negligible level. In particular, we look at all possible 
photon combinations of the four photons and reject any event in which the invariant 
77 mass for both pairs of photons reconstruct near the tt° mass. 

Kl — > vr°7r 7r° decays constitute the largest source of background to Kl — > 7r°77. 
These events can contribute to the background in three distinct ways. First, two 
photons from the Kl — ► 7r 7r°7r decay can miss the calorimeter, with the remaining 
four photons interacting in the Csl calorimeter. Another possibility is that one photon 
misses the calorimeter while two photons overlap in the calorimeter. Finally, all 
six photons can hit the calorimeter, with four of the six photons from the Kl —> 
7T 7r 7r° decay overlapping in the calorimeter, producing four separate clusters in the 
calorimeter. 

The first source of Kl —> 7r 7r°7r° background can be reduced by first eliminating 
events with signals in the photon vetoes. In addition, one can improve the signal 
to background ratio by cutting on the reconstructed z position of the event. When 
the energy in the calorimeter is less than the kaon energy, the event reconstructs 
downstream of its true decay position. This can be seen in Figure El where the signal 
events are relatively flat, while the background events show a large enhancement at 
the downstream end of the detector. 

To get rid of events in which photons overlap in the Csl calorimeter, we define a 
photon shape variable. This variable uses the 3x3 array of Csl crystals containing the 
core of the shower and compares the energy distribution to an ideal energy distribution 
determined from Monte Carlo. This shape variable can be seen in Figure El The top 
plot shows the photon shape variable for Kl —>■ vr 7r° events for both the data and the 
Monte Carlo, and shows good agreement between the two. The bottom plot shows 
the data with the Kl — ► 7r°7r 7r° Monte Carlo overlaid. As can be seen, the signal 
events tend to peak at low values of the shape variable, while the 37r° background is 
relatively flat. We require that the shape variable be less than 1.8, which removes 
the majority of the 3n° events. 

Because of the discrepancy between the NA48 and KTeV published results on 
Kl — > 7r°77, we reexamined the data used in our previous result. [TU] We found that 
we underestimated the impact of a disagreement between the data and the Monte 
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Carlo photon shape variable. This led to an underestimate of the background in 
our final sample. With an improved simulation of the photon shape variable, our 
background estimate nearly doubles, and the Kl — > vr°77 branching ratio decreases. 

After all cuts, the background is dominated by Kl — > 7r 7r°7r° events, with the 
total background constituting approximately 30% of the final events. To check that 
we understand the background level, we examine the data/MC overlays of the m n o 
distributions. Any underestimate of the 37r° background would manifest itself as 
a data/MC mismatch in the tails of the m n o distributions. As shown in Figure EJ 
at small values of the photon shape variable, the tails agree quite well in the m n o 
distributions. 

The systematic errors associated with this measurement are shown in Table H3 
The main sources of uncertainty stem from understanding the 3ir° background and 
its normalization. Other sources of systematic uncertainty come from our knowledge 
of the photon veto system, the acceptance determination and external factors such 
as the measured Kl —> 7r°7r° branching ratios. The total systematic uncertainty is 
2.9%. 



Source of Uncertainty 


Uncertainty (%) 




ay dependence 


1.5 




37T° background 


1.3 




MC statistics 




1.0 


Normalization 




0.9 


Photon Shape 


1.1 




Tracking Chambers 


0.9 




27T° branching ratio 


0.9 




Photon vetoes 


0.9 




Kaon Energy 


0.7 




Decay Vertex 


0.4 




Total 


2.9 





Table 1: Systematic errors for the Kl — > 7r°77 measurement 

After all cuts have been implemented, we find 1982 events with a background of 
601 events. The reconstructed 77 mass is shown in Figure El The distinctive 77 
shape results from coupling of the 77 system to two virtual pions, and peaks around 
320 MeV/c 2 . We determine the K L — > 7r°77 branching ratio to be: BK(Kl — > 7r°77) 
= (1.30 ± 0.03 ± 0.04) x 10~ 6 , where the first error is statistical and the second error 
systematic. As shown in Figure El this result is compatible and competitive with the 
result from NA48. Our new result supercedes the previous KTeV measurement of 
BK(Kl — > 7r°77), which was in nominal disagreement with the NA48 result. 
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4 Measurement of Kl —> tt e + e 7 

The decay Kl —> 7r°e + e _ 7 produces two charged tracks in the spectrometer and three 
photons in the Csl calorimeter. The three photons can be combined in three different 
ways to form a neutral pion. We choose the combination that has the best 7r° mass. 
The neutral vertex is used to determine the decay position rather than the charged 
vertex due to its better resolution. 

Like the K L — > 7r°77 decay, the main backgrounds to this decay are from K L — > 
7r°7r° and Kl — > 7r°7r 7r°. The difference is that one of the tt° undergoes Dalitz decay 
to e + e~7. To help reduce the background from Kl — > 7r 7r° decays, we formed a 
neural net variable using the reconstructed masses m 77 and m e + e - 1 from the second 
and third best combinations. We define the second and third best combinations by 
how far from the nominal ir° mass the 77 combination reconstructs. The neural net 
variable can be seen in FigurelHl In this plot the Kl —> tt°tt events are well-separated 
from the signal Kl — ► 7r°e + e~7 events. We require that the neural net variable be 
greater than 0.5 

To reduce K L — > 7r 7r 7r° backgrounds we also require that the photon shape 
variable, defined above, be small for each photon candidate. We can use the decay 
kinematics to help reduce the background from Kl — ► decays. We define a 

variable p 2 L , which is the longitudinal momentum squared of the missing 7r° in the 
kaon rest frame. We perform a two-dimensional cut on the p L variable versus the 
three photon invariant mass, m 7Tr This distribution is shown in Figure |H1 There 
is good separation between the signal events, shown in blue, and the Kl —* 7r 7r°7r° 
background. We use a polynomial to define the cut shown in Figure 

After applying all cuts, we find 139 events over a background of 14.4 events as 
shown in Figure We reconstructed 80,445 Kl —> 7r°7r° events used for normal- 
ization. This allows us to determine the branching ratio for Kl — > 7r°e + e _ 7 to be 
BR(K L -> vr°e + e-7) = (1.90 ± 0.16 ± 0.12) x 10~ 8 where the first error is statistical 
and the second is systematic. The results from the 1997 and 1999 data sets are shown 
separately in Figured along with the previously published result on Kl —> 7r°e + e~7. 
The published result and our new 1997 result use the same data set. Although the 
branching ratio results for these two results differ, the two results are statistically 
consistent with each other. We expect the two results to differ because of differences 
in selection criteria and calibration constants. As in the Kl — ► 7i~°77, the normaliza- 
tion mode branching ratio has decreased by approximately 8% from the value used in 
the published KTeV result. The systematic uncertainties are listed in Tabled The 
largest systematics come from the limited Monte Carlo statistics and the dependence 
of the result on the parameter ay. 
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Systematic Error (%) 

MC Statistics 4~2 

ay dependence 3.8 

K L and tt° BR 2.8 

3tt° bkg 0.8 

acceptance 0.4 

27T° background 0.1 

Total OA 



Table 2: Systematic errors for the Kl — > ir°e + e 7 measurement 



5 Search for Kl —> n n 7 

In this decay we use the 7r° Dalitz decay for one of the neutral pions. We chose this 
because the fully neutral mode trigger was heavily prescaled, whereas the KTeV two- 
track trigger was not. The final state consists of an e + e~ pair, and three photons in 
the Csl calorimeter. As in the other analyses discussed in this paper, there are three 
possible combinations of photons for each event. Again, we choose the combination 
that reconstructs closest to the 7r° mass. 

The main background to this decay comes from Kl —> 7r 7r°7r° events, since Kl — ► 
7r°7r° cannot contribute to the background without the addition of accidental particles. 
Two variables are effective in reducing the background to this mode. The first is to 
employ the photon shape variable described above. This reduces events in which two 
of the photons overlap in the calorimeter. The second varible used is the missing 
momentum in the kaon rest frame, p\. Cutting on this variable significantly reduces 
the Kl — > 7i 7i°n background. 

A unique feature of this analysis is that the normalization mode is not fully re- 
constructed. Rather we use Kl — > 7r°7r 7r events in which one of the photons passes 
through one of the beam holes in the Csl calorimeter. This photon is kinematically 
constrained using the known kaon mass. 

The final candidates are shown in Figures |Hl where the reconstructed kaon mass 
is shown along the x-axis and the kaon transverse momentum squared, Pj>, is shown 
on the i/- axis. For the 1997 sample we chose a square signal region, while for the 1999 
sample we formed a likelihood from the product of the and p\ probabilities. After 
all cuts we found no events in the 1997 sample and one event in the 1999 sample. The 
probability of this event to come from background is determined to be approximately 
10%. 

Using our Kl — > tt 7t°tt° normalization sample, we determine an upper limit of 
BK(Kl — > 7r°7r°7) < 2.32 x 10~ 7 . This constitutes an improvement of about 20 over 
the previous limit from NA31. However, this limit is well-above the predictions for 
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this decay. 



6 Summary and Conclusions 

The KTeV experiment has presented three new results on Kl — > 7r°77, Kl — > 
7r°e + e~7 and Kl — > 7r°7r°7. The first result is competitive with the world's best 
result from NA48, while the other two represent the world's best measurements on 
these two decays. The measured branching ratios are both inconsistent with the 
0{p A ) predictions from chiral perturbation theory, and consistent with 0{p & ) chiral 
perturbation theory. 
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Figure 3: The photon shape variable for — > 7r°7r° (top) and —> 7r°77 candi- 
dates (bottom). In the top plot, the dots are the data and the red histogram is the 
Monte Carlo. In the bottom plot, the blue histogram represents the Kl — > 7r 7r 7r 
background. 
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Figure 4: The reconstructed 7r° mass distribution for candidate events. The dots are 
the data, the red histogram the sum of the Monte Carlo, and the blue histogram the 
background Monte Carlo. The three plots represent different regions 
of the photon shape variable. 
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Figure 5: The top plots show the final 77 mass distribution for Kl — > tc 77 candi- 
dates before and after background subtraction. The bottom plot shows the branching 
ratio results for the previous KTeV result, the NA48 result and our new KTeV result. 
The results have been rescaled using the latest PDG value for Kl — > 7r°7r°. 
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Figure 6: The variables used to reduce the Kl — > 7r°7r° (top) and Kl 
(bottom) backgrounds. The top plot shows the neural net variable for the Kl — ► 7r°7r° 
background (red) and the Kl — > 7r°e + e _ 7 signal Monte Carlo. The bottom plot 
shows the missing momentum versus the 777 mass distribution for Kl — > 7r 7r°7r° 
background (red) and signal Monte Carlo (blue). 
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Figure 7: The e + e~777 mass (top) for the combined 1997 and 1999 data sets. The 
dots are the data, the blue histogram is the sum of the signal plus background Monte 
Carlo, and the yellow histogram represents the sum of the Kl — > 7r°7r° and Kl — ► 
7r 7r 7r backgrounds. The bottom plot shows the branching ratio results for the KTeV 
PRL (red) and our new 1997 (blue) and 1999 (green) analyses. The results have been 
rescaled using the latest PDG value for Kl — > 7r°7r°. 
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Figure 8: The p\ versus kaon mass distributions for the 1997 (top) and 1999 (bottom) 
samples after all selection cuts have been applied. The red line represents the signal 
region. 
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